Quantum flutter of supersonic particles in one-dimensional quantum liquids 
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The non-equilibrium dynamics of strongly correlated many-body systems exhibits some of the 
most puzzling phenomena and challenging problems in condensed matter physics. Here we report 
on essentially exact results on the time evolution of an impurity injected at a finite velocity into 
a one-dimensional quantum liquid. We provide the first quantitative study of the formation of the 
correlation hole around a particle in a strongly coupled many-body quantum system, and find that 
the resulting correlated state does not come to a complete stop but reaches a steady state which 
propagates at a finite velocity. We also uncover a novel physical phenomenon when the impurity 
is injected at supersonic velocities: the correlation hole undergoes long-lived coherent oscillations 
around the impurity, an effect we call quantum flutter. We provide a detailed understanding and an 
intuitive physical picture of these intriguing discoveries, and propose an experimental setup where 
this physics can be realized and probed directly. 
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Quantum environments are known to be capable of 
drastically altering the properties of embedded parti- 
cles, notable examples being the formation of polarons in 
solid-state systems [I] , Kondo singlets in systems with lo- 
calized impurities [2], and quasiparticles in Fermi liquids 
[3J . The study of these phenomena has typically been car- 
ried out assuming the dressing of the particle is in equilib- 
rium [3] , however recent experiments have begun address- 
ing nonequilibrium phenomena associated with the for- 
mation of these strongly correlated states [MZ] • The the- 
oretical analysis poses one of the most formidable chal- 
lenges in condensed matter physics, requiring accurately 
capturing the dynamics of strongly interacting quantum 
phases of matter [3J H2 110) . In this work we provide an 
essentially exact numerical study of the formation of a 
correlation hole around an impurity injected into a one- 
dimensional gas of hardcore bosons, also known as the 
Tonks- Girardeau (TG) gas [HJH2]. Intriguingly, we find 
the most striking features when the particle is injected 
supersonically. The physics of fast particles is responsible 
for a rich variety of phenomena, such as flutter in aerody- 
namics, Cerenkov radiation |13j . and brcmsstrahlung [14] . 
Our results provide an example of new physics induced 
by supersonic motion in a non-relativistic quantum sys- 
tem. Previous works on fast propagation in Bose gases 
assumed either a weakly coupled gas described by a set 
of noninteracting Bogoliubov excitations [15l - f2"0] . or a 
strongly interacting system treated within a low-energy 
effective field theory approach [2TM23J . In this paper we 
find that there are novel features that require both the 
strong coupling regime and a high energy impurity, thus 
going beyond the regime addressed in previous works. 

Our main observations in tracking the fate of the im- 
purity injected into a ID quantum liquid are twofold. 
Firstly, the injected particle forms a strongly correlated 
state with the quantum liquid that does not come to a full 



stop, instead it reaches a steady state which propagates 
at a reduced velocity. We provide a study of the for- 
mation of these dissipationless propagating states, and a 
protocol for their generation by direct particle injection. 
Secondly, if the impurity is injected at a supersonic veloc- 
ity, the correlation hole around the impurity undergoes 
pronounced oscillations. We call this phenomenon quan- 
tum flutter in analogy with supersonic flutter in aerody- 
namics, which also arises from nonlinear interactions of a 
fast object with the background medium. This quantum 
flutter is due to the formation of an entangled many- 
body state, whose coherence is long-lived. Recent work 
has shown that strongly coupling a particle to a bath 
can lead to non-Markovian dynamics and the possibility 
of coherence surviving for long times [231 HI] ■ Quantum 
flutter provides an example of a quantum system taken 
far out of equilibrium whose relaxation shows striking 
quantum coherent effects that go beyond a hydrodynam- 
ical description. We propose a direct experimental re- 
alization of this physics in a cold atomic setting, where 
impurity physics in Tonks-Girardeau gases has already 
been realized [TJ [26l [27]- But first, we develop a detailed 
understanding of the intricate and intriguing physics of 
the formation of the dissipationless current-carrying state 
and subsequent quantum flutter. 



Physical system and correlation hole formation 

Our system consists of an impurity with flipped spin of 
mass 7Ti j, interacting via a contact interaction with a one- 
dimensional TG gas of hardcore bosons, of mass m\. At 
present, we assume the masses are equal: m^ = m^ = m. 
In the absence of the impurity, Girardeau [12] showed 
that the TG gas can be mapped to a fully polarized non- 
interacting Fermi gas, thus solving an interacting many- 
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FIG. 1. Correlation hole formation and quantum nutter: Formation of the correlation hole and emission of the wave 
packet in the background gas, for interaction strength 7 = 5 and injected momentum Q — 1.35fc.F- a, Time evolution of the 
density distribution around the impurity, L{pf(x,t)p±(0,t))/p, for an impurity with initial momentum Q = 1.35fcF pointing 
in the x > direction. From t — to about t = 5£f, the correlation hole forms. Simultaneously, a large wave packet in the 
background gas forms (red ridge). Friedel-like oscillations in the spatial direction are visible outside of the correlation hole. 
We highlighted the line at t = 25£f outside the correlation hole in red, where the oscillations are clearly visible. In blue we 
highlighted two lines at constant positions (xp = —0.85 and xp = 0.93) on the peaks bordering the correlation hole, along 
which oscillations, which we call quantum flutter, are visible, b, Front view of the correlation hole. Here we see oscillations 
inside the correlation hole, which correspond to the quantum flutter. We highlighted lines of constant height {xp = 0.6 and 
0.8), which are oscillating back and forth in space as a function time, thus the correlation hole is 'fluttering', c, Change in 

momentum distribution in the background gas after the formation of the correlation hole: L\ (hkf(t = 5£f)} — {fik-[(t = 0)} ) 

where hut = cjLcfcf. A narrow peak in the momentum distribution of the majority particles is formed (indicated by the red 
arrow), corresponding to the emitted wave packet. Also, we see a depletion of {hhf) for \k\ < kp, from the correlation hole 
being delocalized across the Fermi sea. Note that (hkf) is defined in the fermionic language. 



body quantum problem exactly. We have now extended 
this argument to the problem with an impurity: because 
the impurity is distinguishable from the TG bosons, the 
TG bosons can be fermionized independently of the im- 
purity. 

Since we are describing results for an impurity in a 
TG gas and a fully polarized Fermi sea simultaneously, 
we can adopt the notation from the fermionic system. All 
the 'fermionic' quantities we describe have direct analogs 
in the TG gas. For example, we define the 'Fermi mo- 
mentum' kp of the TG gas (or the Fermi sea) through 
its density p: kp — ^P- We will refer to either the TG 
bosons or the fullly polarized fermions as the background 
particles, and give them all spin up. The impurity is then 
a down spin. We define fermionic creation (annihilation) 
operators: c k {eta) creates (annihilates) a fermion of 



momentum k and spin a (a =\ or J,). The ground state 
of the background gas becomes simply a fully polarized 



Fermi sea: |FS) = IlifcKfc c Itl^' wuere FS stands for 



Fermi sea and |0) is the vacuum. 

We call xi the position of the impurity, and x%, . . . , xn 
the positions of the N particles in the background gas 
(where N is odd). The Hamiltonian of the system is 
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where P^ (Pjf) is the momentum and m^ (m^) the mass 
of the impurity (i — th background particle); g is the 
coupling between the impurity and the bosons; p a (x) — 
ipl(x)tlj l j(x) is the density operator of spin a, where tp\{x) 
(if} a (x)) creates (annihilates) a fermion of spin a at po- 
sition x. As stated previously we set m^ = m-f = m for 



now. We parametrize the interaction strength between 
the impurity and the background particles with the di- 
mensionless interaction parameter 7 = mg/p. 

To obtain the physics we are describing, we start the 
impurity in a plane wave, thus our initial state is 



i*> = 4,,i F s> 



(2) 



where Q is the initial impurity momentum. We let this 
state evolve in time under the Hamiltonian in Eq. [I] and 
then calculate the expectation value of the quantities we 
are interested in. 

Because this system is integrable, obtaining the many- 
body eigenstates of the problem reduces to solving the 
Bethe-Ansatz equations relevant to this problem 153] . 
Therefore, after obtaining all the eigenstates, one can 
calculate the expectation value of any operator, by in- 
troducing a complete set of eigenstates, and using mas- 
sively parallelized computing resources (see Supplemen- 
tary Information Sections S.l through S.4 for details). 
To probe the dynamics in our system, we calculate the 
density distribution of the quantum liquid around the im- 
purity. This corresponds to calculating the expectation 
value j- (pi(0,t)pf (x, t))L/vF, which measures the den- 
sity of the quantum liquid a distance x away from the 
impurity (we multiply by L and divide by the Fermi ve- 
locity Vf = hkp/m such that this quantity is 1 at time 
t = 0). We plot this quantity in Fig. [Ik for coupling 7 = 5 
and initial momentum Q = 1.35kp. We define a Fermi 
time as the inverse Fermi energy: tp = 1/Ep = 2m/ k p. 
We see that for t from to 5tF, the correlation hole is 
forming, and a narrow wave packet in the background 
gas is being emitted. This last statement can be ver- 
ified by calculating the momentum distribution in the 
background gas, (n^(i)). We plot the change in this 
quantity in Fig. [it at t = 5tp, and indeed find a narrow 
wave packet peaked at some momentum close to Q (Note 
that (rifc-|-) only corresponds to the physical momentum 
distribution in the fermionic case |ll|h 



Phenomenology of the quantum flutter 

In the pictures of the density distribution of the quan- 
tum liquid around the impurity, we saw hints of the quan- 
tum flutter as evidenced by oscillations in time inside the 
correlation hole (see Fig. [lb). To see these oscillations 
in more detail, we plot the time derivative of the density 
distribution j;Ji(A|.(0, t)pf(x,t))L/vF inside the correla- 
tion hole (meaning within an interparticle distance 1/p 
from the impurity). We see that the correlation hole is 
moving back and forth around the impurity. We can find 
a pronounced signature of this effect in another quan- 
tity, the time dependence of the momentum of the im- 
purity (Pj,(t)) which we plot in Fig. [2b. The expected 
impurity momentum initially decays quickly during the 



formation of the correlation hole, but the decay even- 
tually abates, leaving the impurity with a final momen- 
tum equal to a sizeable fraction of the Fermi momentum. 
After this initial drop, the impurity momentum shows 
pronounced oscillations at a fixed frequency. The first 
phenomenon, saturation of momentum loss, happens for 
any nonzero initial momentum, but the oscillations only 
become pronounced when Q is close to or larger than kp. 
The oscillations in the impurity momentum and the oscil- 
lations of the correlation hole are in phase, as evidenced 
in the time derivative of the first moment of the density 
distribution of the quantum liquid around the impurity: 

L J^ ! dxx(pi(0,t)p^(x,t))/vp, see Fig. 2;. This quan- 
tity corresponds to the classical momentum of a density 
distribution of the form (/5^(0, t)p^(x, £)), and it oscillates 
approximately in phase with the oscillations of the impu- 
rity momentum. This result suggests that the oscillations 
are due to a momentum exchange between the impurity 
and its correlation hole. 



Physical mechanism behind the quantum nutter 

Now that we have a complete picture of the different 
components in the time evolution, we detail the physical 
mechanism underlying the quantum flutter: the system is 
undergoing coherent oscillations between two families of 
states which we call 'exciton-like'and 'polaron-like' (see 
Fig. [3h) . Therefore the fluttering corresponds to a quan- 
tum beating between two families of states. That this 
quantum superposition is coherent for such a long time 
is quite remarkable, and is related to the fact that the po- 
laron and exciton dispersions are relatively flat at strong 
coupling. 

To provide insight into the nature of the dynamical 
processes and motivate the intuitive physical picture de- 
scribed above we supplemented the Bethe Ansatz analy- 
sis with a variational approach based on a restricted set of 
wave functions which have been used succesfully in previ- 
ous works to study the ground state properties of impu- 
rities in cold atomic systems 4, 30, ST], and are capable 
of capturing the states that we conjecture are respon- 
sible for quantum fluttering. That this approach agrees 
quantitatively with the Bethe Ansatz results (see Supple- 
mentary Information Section S.5) strongly supports the 
qualitative picture we will now describe. We emphasize, 
however, that all quantitative results presented in this 
work (except for the mass imbalance results in Fig. p| 
were obtained from the full Bethe Ansatz equations. 

Our intuitive picture behind quantum flutter relies on 
two types of states, an exciton and a polaron, which we 
describe qualitatively using variational wave functions. 
The exciton state with total momentum K is the lowest 
energy state composed of a Fermi sea, an impurity and 
a hole. In its simplest incarnation it has the following 
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FIG. 2. Properties of quantum nutter: Properties of the quantum flutter, with parameters 7 = 5 and Q = 1.35/cf. We 
observe oscillations of the correlation hole around the impurity, and of the momentum of the impurity, which are in phase 
and consistent with the picture of the impurity exchanging momentum with its correlation cloud. This suggests that coherent 
quantum mechanical processes are leading to a collective mode of the system, which we call quantum flutter, a, Plot of the 
time derivative of the density distribution of the background gas within the correlation hole: jr-^{p±(0,t)p^(x,t))L/vF with 
x within an interparticle distance p _1 of the impurity. Pronounced oscillations of the correlation hole are observed at a fixed 
frequency, b, Time evolution of the expected momentum of the impurity. We notice two main features: first the momentum 
does not decay to zero. Second, the momentum shows oscillations at a fixed frequency around the saturation value. These 
oscillations are in phase with the oscillations of the correlation hole, c, Comparison of the time derivative of the first moment 

of the density of the background gas inside the correlation cloud, L-^ J p _ 1 dxxpi(0,t)p^(x,t)}/vF, and the oscillations of the 
impurity momentum. If the background gas were classical, this time derivative of its first moment would correspond to its 
momentum. The oscillations of the two quantities are almost in phase and of the same order of magnitude, substantiating the 
claim that the quantum flutter corresponds to the impurity exchanging momentum with its correlation hole. 



form: 



\Exc(K)) = Q2. 
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We call E(Exc(K)) the energy of the exciton at momen- 
tum K . We remind the reader that Eq. [3]gives a descrip- 
tion of the exciton in a truncated Hilbert space. A full 
description of this state, obtained from the Bethe-Ansatz 
equations, will contain an infinite number of particle- 
hole pair excitations. However, this level of description is 
enough to obtain quantitative agreement with the Bethe 
Ansatz calculations, as we describe later. 

A polaron state is a dressed impurity state, i.e. it 
consists of an impurity on top of a Fermi sea, dressed 
by particle-hole excitations of the Fermi sea. We call 
\Pol{K)) the ground state of a system made up of a Fermi 



sea and an impurity with total momentum K, which, 
allowing at most one particle-hole pair excitation, is given 
by [30] 



\Pol(K)) = (0 



ff) 4++E7ff4- fe+g +4t c «t)l FS >- 

kq 



(4) 

We call E{Pol(K)) the energy of the polaron at momen- 
tum K. As for the exciton, in the full Bethe Ansatz 
solution the polaron will contain an infinite number of 
particle-hole pair excitations, but capturing the polaron 
at this level already gives quantitatively accurate results. 
Now, once the system has emitted a wave packet in the 
background gas, it leaves behind a hole that the impurity 
can interact with. Therefore the impurity and hole can 
form an exciton, and the wave function of the system 
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FIG. 3. Physical picture behind quantum nutter: An 

intuitive picture behind the quantum flutter involves the sys- 
tem forming a coherent quantum superposition of two fam- 
ilies of states. The energy difference between the van Hove 
singularities, where the density of states diverges, in these 
families of states is exactly equal to the oscillation frequency 
of the quantum nutter, which becomes a quantum beating 
effect in this picture, a, The families of states behind the 
quantum nutter, as discussed in the text. The impurity 
can bind to the hole left behind by the emitted particle, 
creating an exciton. Alternatively, it can also not bind to 
the hole and formed a state dressed with particle-hole pairs, 
called a polaron. Therefore, the two families of states are 
c k AExc(Q — k)) (Exciton+Particle) and c k *.c q f\Pol(Q — k+q)} 
(Polaron+Hole+Particle). b, Stars: oscillation frequency 
hw osc = h2n/r osc in the impurity momentum where t osc is 
the period of the oscillations seen in Fig. [2b. Green line: 
energy difference between CqJExc(0)) and Cq^Coi\PoI(0)} , 
which equals E(Pol(0)) — E(Exc(0)). They agree to within 
1%. 



would be 



Exciton+Particle state: c k AExc{Q — k)) (5) 

where k is the momentum of the emitted particle. An- 
other option is that the impurity does not bind with the 
hole, instead it forms a polaron with the Fermi sea. We 
then have a polaron, a hole and a particle in the back- 
ground, giving a wave function of the form 

Polaron+Particle+Hole state : c k *.c q ^\Pol(Q — k + q)) 

(6) 
See Fig. [3£i for a schematic representation of these two 
possibilities. 

One can interpret the quantum flutter as being due to 



the system being in a superposition of these two possibili- 
ties: Exciton+Particle and Polaron+Particle+Hole. The 
particle has momentum close to the total momentum Q, 
and the hole is a deep hole, meaning that its momen- 
tum is close to zero. Consider now the case where the 
particle has momentum exactly equal to Q, and the hole 
momentum exactly equal to zero. Then our two states 
are c*qAExc(Q)) and Cq+cq^\PoI(0)). Now, we can cal- 
culate the energy of \Exc(0)) and \Pol(0)) exactly using 
Bethe Ansatz (this means that Eqns. [3] and El are now 
dressed with an arbitrary number of particle-hole pair ex- 
citations). Using this quantity, we can look at the differ- 
ence in energy between Cq^co^\PoI(Q)) and CqJExc(0)) , 
which is E(Pol(0))—E(Exc(0)), because both states have 
a Fermi sea and a particle at Q so those energies cancel 
out, and the hole at zero momentum carries no energy. In 
Fig. [3b, we plot this difference in energy and compare it 
to the frequency of the quantum flutter hu osc = h2ir/T osc 
where t osc is the period of the oscillations in impurity 
momentum. The two quantities match to within 1%, 
which is consistent with our physical picture of quantum 
flutter arising from quantum beating. Note that in cal- 
culating the energy we used the fact that the emitted 
wave packet is spatially separated from the impurity and 
therefore does not interact with it. Also, we assumed 
the deep hole weakly interacts with the polaron, as the 
polaron is mostly dressed close to the Fermi points. 

As we mentioned earlier, the emitted particle does not 
have momentum exactly equal to Q, and the hole on top 
of the polaron is not exactly at zero momentum, instead 
they have some spread. However, at stronger coupling 
the polaron and exciton have increasingly flat disper- 
sions. Therefore as long as the hole on top of the polaron 
is close to zero momentum, then the energy difference 
between c k ^c^\Pol(Q — k)) (setting the hole momentum 

to zero for now) and c k JExc(Q — k)) will depend only 
weakly on k. Thus the emitted particle can form a wave 
packet without disturbing the oscillation frequency of the 
flutter much. However, forming a superposition of states 
with different energies will lead to damping of the os- 
cillations, which is indeed what we observe. From this 
picture we can also understanding why the oscillations 
only appear when Q is of the order of or larger than kp: 
Q = kp is the minimum momentum necessary to be able 
to create a particle in the background gas (whose mo- 
mentum has to be larger than kp) and an exciton at zero 
momentum. That the oscillations become so pronounced 
is then due to the van Hove singularity at the top of the 
exciton branch (i.e. in E(Exc(K)) at K = 0). 

As mentioned earlier, the picture of quantum flutter 
given above is qualitative, but in fact a variational ap- 
proach that can capture the families of states in Eq. [3] 
and Eq. pleads to results which are in quantitative agree- 
ment with the results from Bethe Ansatz (see Supplemen- 
tary Information section S.5 for details). This strongly 



supports the physical understanding behind quantum 
flutter that we expound, of a long-lived coherent super- 
position of exciton-like and polaron-like states. 
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FIG. 4. Dependence of oscillations on initial mo- 
mentum: Average impurity momentum (Pj_(t)} as a func- 
tion of time, for 7 = 5, for different initial momenta: (red) 
Q — 1.05kp, (orange) Q — 1.35/s.p, (green) Q — 1.7kp, (blue) 
Q — 2kF- Inset: zoom in on the oscillations. The structure of 
the oscillations depends only weakly on the initial impurity 
momentum. This implies (see text) that if the impurity was 
in a wave packet state peaked at high enough momentum, the 
average momentum of the impurity would flutter. 



Experimental consequences 

We propose several experiments to detect the physics 
described above. For the most direct realization of it, 
consider a neutral bosonic atom that is confined to a 
quasi- ID geometry and made to have strong repulsive 
interactions via a Feshbach resonance, thus realizing a 
simulation of a Tonks-Girardeau gas [7, 26, 27, 33J . This 
atom must have at least three internal hyperfine states 
|1), |2) and |3), such that |1) and |3) strongly interact, 
while the pairs |1) and |2), and the pairs |2) and |3) 
weakly interact. Start the system with most atoms in 
the |1) state, and a few in the |2) state. Now use a two- 
photon Raman pulse to excite the |2) atoms into a |3) 
state with a given momentum. Then use time-of-flight 
measurements to map out the time dependence of the 
momentum distribution of the |3) particles. The aver- 
aged momentum should see the quantum flutter. 

Typically impurity experiments use localized RF 
pulses to create impurity wave packets. From our re- 
sults for (Pj.(t)), one can immediately obtain the results 
for (Pj_(i)) for an impurity initially in a wave packet, i.e. 
an initial state of the form 2fc a fc c L|FS) where at is 
the Fourier transform of the wave packet. Namely the 
resulting (Pj,(i)) can immediately be obtained from our 
calculations by averaging the results obtained for plane 



The reason is that the Hamiltonian conserves total mo- 
mentum, so the different momentum components of the 
wave packet evolve independently of each other. In Fig. [4] 
we show that (GS\ck^P^(t)c\, \GS) is only weakly depen- 
dent on k after the initial decay if k > kp, therefore the 
quantum flutter would not be washed out if the wave 
packet was peaked around high enough a momentum. 
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FIG. 5. Dependence of oscillations on mass imbal- 
ance: Time dependence of the impurity momentum (Pj.(i)) 
for initial momentum Q = 1.05kp, interaction strength 7 = 5, 
and different mass imbalances r — m^/rrif, from top to bot- 
tom: r = 0.8,1,1.2,1.4. These results are obtained from 
the variational approach discussed in the text. The r = 1 
case agrees quantitatively with the result obtained from the 
Bethe Ansatz approach (see Supplementary Information Sec- 
tion S.5). We see that the saturation of momentum loss and 
quantum flutter exist away from the integrability point. For 
r < 1 the quantum flutter gets strongly damped, while for 
r > 1 the damping depends only weakly on r. 

Whether quantum flutter is present in other models is 
an important question, which we have begun answering. 
As discussed above, we have a variational approach that 
agrees with the Bethe Ansatz results, and therefore al- 
lows us to explore non-integrable models, ft is known 
that the equilibration of integrable and non-integrable 
models is qualitatively different, due to the infinite set 
of conserved quantities in the integrable case pH] . Non- 
integrable models, on the other hand, are expected to 
relax to a local equilibrium described by hydrodynamics. 
In our case, we varied the ratio r — m^/m-f of the mass 
of the impurity and of the background particles, making 
the system non-integrable. As shown in Fig. [5] we find 
that the saturation of momentum loss and subsequent 
quantum flutter are still present when the mass ratio is 
changed. For r < 1 the quantum flutter becomes more 
strongly damped, while for r > 1 the damping is weakly 
dependent on mass ratio. 



Discussion 



We have provided an example of new physics arising 
from the injection of a supersonic particle into a many- 
body quantum system. In many physical and biologi- 
cal systems |35l I36| an important question is whether a 
correlated quantum state can travel through an environ- 
ment while maintaining its coherence. The phenomenon 
of quantum flutter provides an example of the forma- 
tion of an entangled state that is propagating through 
the many-body system and remaining coherent for long 
times. Furthermore, we showed that this dynamically 
generated and protected coherent state is not a singular 
feature of the integrable point, and we find similar fea- 
tures as we perturb away from integrability. The most 
important question that remains from this work is which 
physical systems exhibit effects similar to quantum flut- 
ter, and more generally whether quantum systems taken 
far out of equilibrium exhibit quantum coherent effects, 
which go beyond a hydrodynamic picture of local relax- 
ation. We have begun answering these questions, and 
believe that state-of-the-art theoretical and experimental 
methods for studying nonequilibrium quantum dynamics 
should be able to shed more light on this topic. 
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SUPPLEMENTARY INFORMATION 

THE MODEL AND POLARON 

TRANSFORMATION TO THE IMPURITY 

FRAME 

We are considering a one-dimensional Fermi gas of 
atoms which we call spin up, and one impurity with spin 
down (in other words, we fermionized the TG gas). The 
second-quantized representation of the Hamiltonian is 

11 = ~ f L dx J2 id 2 Ji(xM a (x) 



operator S = (S x , S y ,S z ): 



2m 



o 



<*=t,l 



+ 9 dx : p f (x)pi(x) : 
Jo 



(SI) 



where p a (x) = tpl l (x)ip a (x), a =t, -I ■ The Fourier trans- 
form 



to = 4f5>~ i9X 4 



(S2) 



where a =f , J,, q 



L ' 



0,±1,±2, ..., gives the 



momentum-space representation of the fermion fields en- 
tering Eq. (Si]). The fermion anti-commutation relations 
in the momentum space are [c aq ,Cg ,] + = 8 a pS qq '. The 
symbol "::" stands for the normal ordering. 

The Hamiltonian (H1J) commutes with the global spin 



[H, S] = 0, S = 



- / (iiV^(i)(T^^(i), (S3) 



where cr = (<J Xl o- y ,a z ) is the vector composed of the 
three Pauli matrices. The spin-ladder operators S± = 
S x ± iS y convert spin-down particles into spin-up, and 
vice versa. Their momentum-space representation is 



5_=^c| fc c Tfc , S+=Sl. 



(S4) 



Another conserved quantity is the total momentum, 

[H,P}=0, P = P t + P l . (S5) 

The model (S[l]) is Bethe-Ansatz solvable. McGuirc 
\S1\ found its wave functions \^f q ) and spectrum E-q, in 
the sector with one spin-down particle, and an arbitrary 
number N of spin-up particles: 

H\%)=E* q \%), ^1%) = ]%), P\%) = q\%), 

(S6) 
where N± — L dx pi(x). Note that we write N for the 
number of spin-up particles, which we also refer to as the 
host particles, instead of N-f to lighten the notation. We 
assume N is odd. 

Despite the McGuire solution, the calculation of 
the far-from-equilibrium impurity dynamics in the 
model (SfTl) still remains a challenge. Call 

K) = cJ,|FS) (S7) 

the initial state where the Fermi sea for N particles is 
defined as (q m = 2i:m/L and 'vac' is the vacuum) 

(N+l)/2 

i fs >= n c i> ac >- ( sg ) 

m=-(iV-l)/2 

The average impurity momentum P q (t) is then 

P q (t) = (m q \P±(t)\m q ), (S9) 

and we refer to it interchangeably as (Pj,(t)) for short- 
hand. It can be represented as 



PS) 



<y U( Eq , q -l:,. :: 



«'(FS|c^|*, 



x{* q \P^ q ){%\cl\FS) (S10) 

using the completeness of the basis of the functions |^g). 
In order to get P q (t) for sufficiently broad ranges of q, 
t, and particle number reliably one should include a vast 
number of terms in the sum on the right hand side of 



Eq. (£10 1. Therefore it is crucial to find a computation- 



efficient representation of the matrix elements entering 



the decomposition (£ 10 1 



To perform this task we find it convenient to use the 
mobile impurity reference frame. The transformation to 
this frame (often called polaron transformation) for an 
arbitrary operator O reads (see, for example, Ref. \S2\ ) 



0->O a = QOQ- 1 , 



oiPf^i 



(Sll) 



where P-|- is the total momentum operator of the back- 
ground particles, and ah is the impurity coordinate oper- 
ator. This transformation acts on the impurity and host 
momenta as follows 



P^Q — -Pf' P-l-Q — P i ~ *\i 



(S12) 



Applying the transformation (£ 11 ) to the Hamiltonian of 
our model, we get 



He 



, )m <te[^(x)]^t(a;) + (n-A) a +3i5t(0)- 

(S13) 
The conservation of the total momentum in the labora- 
tory frame, Eq. (95J), leads to the conservation of P± in 
the impurity frame, [-ffg,PJ = 0. Therefore the func- 
tions l^q), Eq. (£j6J), can be written as 

\%) = e -1 4jo> ® 1/,), (si4) 

where |0) is the state with no particles, and the states \f q ) 
are orthogonal, normalized and form a complete basis of 
eigenstates of the Hamiltonian Hq(c[): 

H Q (q) = ]T \f q )E fq (f q l (f q \f q ) = S ftfi , (S15) 

h 

where 



H Q (q) = - 



., dx[d 2 J\(x)}^(x) + (q-P t ) 2 +gp t (0) 

(S16) 



is obtained by projecting the Hamiltonian (£ 13 1 onto the 
sector with P± equal to q. We stress that \\fr q ) and |/ g ) 
are in one-to-one correspondence, and E^, = Ef . Sub- 
stituting the expression (£14 1 into Eq. (£10) we get 



P q (t) 



q- E e 

JqiJq 



■« E i«- E fO { FS\f q )(f q \P t \fi)(f q \FS). 

(S17) 



SLATER DETERMINANT REPRESENTATION 
OF THE WAVE FUNCTIONS 

A great advantage of the states \f q ) is that in the co- 
ordinate representation these states are Slater determi- 
nants: 



fg(xi 



Y, 



,x N 



./', 



Wl 



det N (4>j(xi)). 



(S18) 



Here Yf is a normalization constant, 



cj>j(x) 



exp < 1 



2^ 
L 



9W 



(S19) 



where j = 1 , . . . , TV and 

, N+l 

9to = 5 E 



9 (-*<*» \'[T* + St 



where 6 = J2t=i @t and 6 t 
8t are 



5 t = h arctan(az t — c), 



(S20) 

a sin St . The phase shifts 
t=l,...,N+l (S21) 



and 



4A 



9L 1 



(S22) 



The rapidities Zt are solutions to the Bethe equations 

cot z t =az t -c, t = l,...,N + 1 (S23) 

and 

o iV+1 



L 



E 



z* = g. 



(S24) 



We stress that it is Eq. (£ 24 ) which ensures that the total 



momentum of the system is q and couples Eqs. (9231) for 
different t = 1 , . . . , N + 1. We will assume that z t are 
ordered in ascending order: 



Z\ < ■■■ < Zjv+1- 
The energy Ef of the state \f q ) is 



E h-J2 



N+l 

E- 



(S25) 



(S26) 



Among the solutions to the Bethe equations (£23) 



there is a special subset corresponding to c = —00. One 
has for the solutions from this subset 



Zt 



nn t , o t = 0, — 



1 



e n + i 



t = l,...,A+l, 

(S27) 

where n t are arbitrary integers. It is thus natural to 
associate Zt with the momenta of a free-fermion problem. 
Let us apply the operator £_ given by Eq. (£J4J) to the 
wave function of N + 1 spin-up fermions: 



1 



1 ql Vn + i tZ1 



-.S-<* 



-\z N . 



JO), 



(S28) 



where q is defined, as usual, by Eq. (£24). Clearly, the 



functions (£28) are indeed eigenfunctions of the Hamil- 
tonian (£[!]). Note also a determinant representation for 
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f q alternative to the one given by Eq. (£18): 



f q {xi,...,x N ) 



1 



y/(N + l) 



J_ e ^X 1 Z 1 ( , t _J_ e ^X X Z N _J_ e J?;X 1 Z N + 1 



1 ^j^^Af-Sl , , , 1 „-^XmZn _J^g x^iV-ZJV + l 

vL vL vL 

1 ••• 1 1 



(S29) 



For c = oo we get the same set of solutions as for 
c = — oo. Therefore only one of the sets should be used 
to avoid double counting. 



DETERMINANT REPRESENTATION OF THE 
MATRIX ELEMENTS 



Using the identity 



1 

m 



dx\ ■ ■ ■ dxN det n[4>j(xi)] det N[fj(xi)] 



det 



N 



dy ct>j(y)<pi(y) 



(S30) 



valid any functions <f>j and (pj, j = 1, . . . ,N we get the 
matrix elements entering Eq. (£ 17) is the form of deter- 
minants of N x N matrices. 

To ensure the normalization (f q \f q ) — 1 of the wave 
function (£ 18 1 the constant Yf should satisfy 



1 f N+1 ffi \ N+1 

\t=i * / *=i 



), (S31) 



where 9 t and 8 are defined by Eq. (£ 20 ) 



The ground state wave function of N free fermions, 
IFSjv), is 



\FS N ) 



1 



'N\ 



:det 



JV 



1 



2uj 



(S32) 



where the momenta u.,-, j = 1, . . . , N, are (we require N 
to be odd) 



N + 1 



.7 = 1, 



,iV. 



(S33) 



For the overlap of the states (£ 18 1 and (£ 32 ) the identity 
(£30) implies 



(FS N \f q )=Y f dct N X 1 (S34) 

where the entries of an N x N matrix X are 

JV+l 



X 



1 



1 yC gt 



Wj - 2TJ 6^Mj- Z f 



The denominators on the right hand side of this equation 
may vanish in the case of c = — oo only. To avoid dealing 
with the resulting singularities it is convenient to use the 



representation (£29) to calculate the overlap with the 
function (£p52l). This gives 



det at X = 



(-1) [ 



Uj = z P . 



(S36) 



N + 1 ' 
and zero otherwise. Here P denotes a permutation of the 



ordered set (£ 25 ) 



Zl, 



• Zn+i — > Zp 1 . 



■i z p* 



(S37) 



and [P] is the sign of the permutation. 

The determinant representation of the matrix elements 

(fMf'g) * 



(f q \P t \f q ) = Y f Y f ,-^det N (y + \Z) 



(S38) 



A=0 



where N x N matrices y and Z are, respectively, 

pL q N+1 

y\ = / dy^vMy) = K(zi, Zj) - ^ E K M> z ') 
J ° U t=i 



% N+1 



3/ N+1 



e 



ffi« + SE^ 



*=i 



■wj 



(S39) 



and 



Zj =-i / dycj) :j {y)dy4>' l {y) 



N+1 



a/ W+l 



i 1 ' r,,J ' e 2L^-«-v~"^ e 



h'^U) - g £ ^(*U) ^E *i K (*t,*i) 

(S40) 



, JV+l 



ee 



t,*'=i 



where 

X(z', z) = g^^-Z-l e'P'-i), K(z, z) = 1. (S41) 

2«(z' — zj 



CONVERGENCE OF NUMERICS 

Now that we've derived how to calculate matrix ele- 
ments like (f q \fg) and (f q \f q ) and (f q \Pf\f q ), we discuss 
how to generate all the states \f q ). One has to generate 
all ordered sets of N + 1 integers {m, . . . , uat+i}. We 
then look for a value of c such that one can solve Eqns. 



j,l = l,...,N. 
(S35) 



(£J23J) and (£24). Not all ordered sets of integers will 



allow for such a solution due to the quantization condi- 
tion 2 J2i z i — qL, therefore one should search through 
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the space of ordered sets of integers for sets that will al- 
low for a set or rapidities that satisfy the quantization 
condition. 

Since the number of many-body eigenstates is expo- 
nential in the system size, the challenge is to generate 
eigenstates that contribute significantly to the calcula- 
tion of a given quantity. Note that although the idea is 
simple, one requires a clever algorithm to generate the 
eigenstates that contribute significantly to the calcula- 
tion. Indeed, it turns out that a set of states that grows 
polynomially with system size is sufficient to obtain accu- 
rate results, making it tractable on massively parallelized 
computing resources. As was shown in other works \S3\ . 
we found that it is convenient to start with the set of inte- 
gers symmetrically ordered around 0, and consider sets of 
integers obtained from this set by taking a fixed number 
of integers out of that set and into integers outside the 
set. In fact, by taking the sets of integers where either 
0, 1, 2 or 3 integers have been taken out of the symmet- 
rically ordered set, one obtains convergent results, as we 
show in this section. This set of of ordered integers grows 
polynomially with N, thus the complexity of the problem 
has become polynomial in system size and tractable. 



We calculate Eq. (E 17) numerically for a finite number 



of particles, N, and for a finite number of intermediate 
states \f q ) and \f) included in the sum. In the present 
section we demonstrate that our choice of N and of the 
most relevant subset of states \f q ) and |/') makes the re- 
sults presented in the paper indistinguishable from what 
one would have in the thermodynamic limit (that is, for 
N = oo) and all intermediate states included in the sum. 
Let us start with discussing the dependence of P q {t) 
on the choice of intermediate states entering the sum in 
Eq. (£|17|), for a given finite N . We define 



Qs = Y^ K/9l FS >| 2 > s = subset of statcs {/<?}■ ( S42 ) 

s 

The completeness of {f q } implies the sum rule £•{/,} = 1 
(i.e. if all states are summed over), therefore we can use 
the devation of 1 — p s from as a check of the convergence 
of our numerics. Let N s be a number of states in the 
subset s. Clearly, p s depends not only on N s but on which 
states are included in s. We plot in Fig. SI the function 
1 — p s versus N s obtained for the states ordered in a 
descending order of the overlap |(/ q |FS)|. 

We then examine how {P\{t)) depends on the number 
of states N s included in the sum in Eq. (£17). A typical 



result is illustrated in Fig. S2. It reflects two important 
facts: (a) our numerics have converged; (b) the error does 
not increase over time, instead it remains bounded. 

We can now examine the dependence of (P±(t)) on N. 
The convergence of our numerics, discussed in the previ- 
ous paragraphs, allows us to separate the effects of finite 
N s from those of finite N. The latter are clearly seen in 
Fig. S3. We see that finite N effects play a role for t 
above some critical value, which increases with N. At 



this critical time one sees revivals, which correspond to 
the emitted particle going around the system. Below this 
time the results are independent of particle number, and 
therefore valid for the thermodynamic limit N = oo. For 
example, Fig. S2b of our paper is a plot for N = oo. 
Indeed, the saturation of momentum loss and quantum 
flutter are obtained before the revivals. 



THE VARIATIONAL APPROACH 

We consider a variational wave function with two 
particle-hole pairs 



(ao(i)cQ + J2 ^«(*) c Q-fc+ Q 4t c 9t 



kq 

+ £ 7fefc ' w '(i) C t Q _ fe+g 4 t c gt )|FS).(S43) 

kk' qq' 

This set of wave functions can capture the states we 
stated were responsible for quantum flutter, Eqs. 3 and 4 
in the main text. We obtain time-dependent variational 
equations by minimizing the expectation value of idt — H 
where 



* = E 



ka 



k 2 

2m n 



9 V" J 



„t 



t 2^ C k-qt C k'+qi C k'lCkt- (S44) 

k,k' ,q 



We therefore solve the time dependent equations from 
considering the variational minimization with respect to 
the variables a (t), fi kq (t), and Jkk'qq'{t)- 



5(*(t)\(id t -H)\*(t)) = 



(S45) 



The resulting coupled equations can be obtained from the 
equations one gets from studying the ground state \S4\ , 
by replacing the energy E replaced with id t . The re- 
sulting equations contain continuum variables k,k',q,q' 
which are discretized, and the equations are solved use 
Runge-Kutta integration. In Fig. S4 we compare the 
results obtained from the variational approach (which is 
in the thermodynamic limit) and the Bethe Ansatz re- 
sult with 48 particles, for 7 = 5 and Q = 1.05kp. We 
find that the absolute difference between the two plots 
remains below 0.02kp. 
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Figure SI. Saturation of the sum rule as a function of the number of states. Shown is the log-log plot of the sum 
rule 1 — g s versus the number of included states N s . The states in the subset s are ordered in decreasing order of |(/ 9 |FS)|. The 
parameters are: number of particles in the background gas TV = 31, dimensionless coupling 7 = 5, initial momentum q = 1.5&F. 
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Figure S2. Convergence of our numerical calculation of (Pj.(£)) using the Bethe-Ansatz approach. 

a, Plot of (P_i{t)} for the number of included states N s = 3 x 10 4 (8 x 10 4 ), which corresponds to a saturation of the sum rule 
g s — 0.994 (0.998). b, Plot of the relative error for these two curves given on the left panel. The input parameters are the 
same as those used for Fig. SI: number of hosts TV = 31, dimensionless coupling 7 = 5, initial momentum q — 1.5/cf. 
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Figure S3. Dependence of our results on the number of particle N in the background gas. We find that finite size 
effects manifest themselves as a revival which is pushed out to longer times for larger particle numbers. Before the revivals, 
the results are independent of particle number. Therefore, since the saturation of momentum loss and quantum flutter occur 
before the revivals, these results are valid for the thermondynamic limit N = oo. a, Plot of (P±(t)) for 7 = 5, initial momentum 
q = 1.5&F and various N: N = 28, 32, 36. b, Zoomed in plots from the left panel. We notice that the curves lie exactly on top 
of each other before the effects or revivals start to set in. 
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Figure S4. Comparison of the Bethe Ansatz results with the variational approach, a, Plot of {P±(t))/kF obtained 
from the Bethe Ansatz approach(full line) and from the variational approach (dotted line), for 7 = 5 and Q = 1.05/sf. The 
Bethe Ansatz results were obtained for 48 particles, b, The absolute difference between the two curves in a is below 0.02fc_F at 
all times. 



